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One of the most exciting discoveries in strongly correlated systems has been the existence of a 
superconducting dome on heavy fermions close to the quantum critical point where antiferromagnetic 
order disappears. It is hard even for the most skeptical not to admit that the excitations which 
bind the electrons in the Cooper pairs have a magnetic origin. As a system moves away from 
an antiferromagnetic quantum critical point, (AFQCP) the correlation length of the fluctuations 
decreases and the system goes into a local quantum critical regime. The attractive interaction 
mediated by the non-local part of these excitations vanishes and this allows to obtain an upper 
bound to the superconducting dome around an AFQCP. 



I. INTRODUCTION 

The study of heavy fermion materials is an exciting 
area in physics motivating sophisticated experimental 
work and giving rise to many new concepts and ideas^. 
In particular the fact that heavy fermions are close to 
a magnetic quantum critical point^ has brought a new 
range of possibilities to this field, both theoretically and 
experimentally. 

On the course of their investigations, as experimen- 
talists aim to reach closer to the AFQCP at even lower 
temperatures, came out the exciting discovery of a su- 
perconducting dome encircling a putative AFQCP^. The 
region of superconductivity in the phase diagram is re- 
stricted to a close neighborhood of the AFQCP. Even for 
the most skeptical it is hard not to admit that in this 
case superconductivity is due to quantum antiferromag- 
netic fluctuations associated with the QCP^. 

The theory of superconductivity mediated by spin 
fluctuations has progressed very much in the last 
decades mostly due to its relevance for high temperature 
sup er conduct ivity^ii^i^i^iS. In these theories, the para- 
magnon propagator describing critical antiferromagnetic 
fluctuations close to an AFQCP can be written in the 
scaling form^. 
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where xs = Xo/ldl is the staggered susceptibility, ^ = 
y^A/\g\ and r = tq^^ the correlation length and critical 
relaxation time, respectively. The quantity g measures 
the distance (in energy scale) to the AFQCP (at ^ = 0), 
A is the stiffness of the spin fluctuations, tq = 1/A and 
the dynamic exponent z = 2. We are interested here in 
quantum phase transitions which occur in three {d = 3) 
or two dimensions {d = 2). Since the dynamic exponent 
associated with an AFQCP takes the value^^ z = 2, the 
effective dimension associated with the antiferromagnetic 
quantum phase transition is dg// = d-\-2. Then, for 2d 
and 3d systems the effective dimension coincides or is 
above the upper critical dimension = 4, respectively. 
This implies that the transition at T = 0, ^ = is de- 
scribed by Gaussian or mean- field exponents^ ^ with log- 



arithmic corrections in the marginal case {d = 2). Then 
for deff > dc the Gaussian exponents, 7 = 1, for the 
staggered susceptibility and u = 1/2, for the correlation 
length turn out to describe correctly the quantum critical 
behavior of the AFQCP. 

The Gaussian free energy close to the AFQCP can be 
written as {ks = 1)^^^, 
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For temperatures T « Tcoh the free energy is given by. 
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where Qc is a cut-off. The coherence temperature, Tcoh = 
l^r^ = 1^1 is that introduced by Continentino et al.^- 
and marks the entrance of the system in the Fermi liquid 
regime. Sd is the surface of a d-dimensional sphere with 
unit radius. The specific heat C/T = —d'^f/dT'^ in the 
Fermi liquid regime, T « Tcohi is easily obtained. 
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in 2d. In the critical regime qd 1, we get that 7 = 
C/T is constant in 3d and logarithmically divergent in 
2di^. We can also define^^^^^ a local limits qd < 1, in 
which case the specific heat is given by. 
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independent of dimension. This result can be obtained 
directly from Eq. [2] neglecting its q-dependence and re- 
placing ^ N. The propagator associated with these 
local spin fluctuations is given by. 
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with xs = Xo/ldl ^ = ^0^^' It is remarkable that 
in spite of the local character of the fluctuations in this 
regime, the system is still aware of the quantum phase 
transition through the dependence of r and xs on g. 
Indeed, in this regime the fluctuations are local in space 
but correlated along the time directions. The theory 
of this regime can be described as a critical theory in 
deff = z = 2 but with Euclidean dimension d = 0. The 
properties of the system for qc(, < 1 have been described 
in Ref.^^. They can all be expressed in terms of a single 
parameter, the coherence temperature^^. 



II. RELATION TO SUPERCONDUCTIVITY 

The q-dependent propagator given by Eq. [T] provides 
an attractive interaction in the singlet channel among 
quasi-particles in neighboring sites. Thus, antiferromag- 
netic paramagnons can give rise to d-wave pairing with 
dx2_y2 symmetry^, for QcS, ^ 1. What about local spin 
fluctuations? 

The spatial dependence of the dynamic susceptibility 
can be obtained from the Fourier transform of the q and 
CO dependent propagator. 
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In the case of local spin fluctuations, the relevant propa- 
gator is given by Eq. [7J such that. 
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Then, 



^exL{uj = 0) = xs- 



X{r) =XsYl e^'-' = 2^xs [ dqdOq' sin^e^^^-^^ 

The interaction among the quasi-particles according to 
Monthoux et al.^ is given by: 

U{r)=-X\{r)Si.S2 (10) 

= -A (-1) o—-nfF\3 1 sm cos — ^ . 



For small r this yields. 
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If we are dealing with a lattice then we have a delta- 
function at the origin. The potential is repulsive at the 
origin and zero everywhere else in the singlet channel. 

Also, it is interesting to obtain the results in the case 
the local propagator is associated with density fluctua- 
tions. In this case^-, 



U{r) = -X\{r) 



and we obtain, 
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where xs in this case is the compressibility. In the case 
of a lattice the interaction is attractive at the origin and 
zero everywhere else in the singlet channel. 

In the magnetic case, Eqs. [11] show that the on-site 
and nearest neighbor quasi-particle interaction mediated 
by critical local spin fluctuations are repulsive in the sin- 
glet channel and do not lead to Cooper pair formation. 
Then, as the system moves away from the AFQCP, the 
correlation length of the spin fluctuations decreases and 
for qc£, ~ 1 the relevant interactions mediated by these 
fluctuations become repulsive everywhere destroying su- 
perconductivity. On the other hand, local charge fluctu- 
ations can still mediate an attractive local interaction. 

The condition 
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puts an upper limit to the region of the phase diagram 
around the AFQCP where superconductivity mediated 
by spin fluctuations can exist. At zero temperature this 
implies that superconductivity survives up to a critical 
coupling \g\s =Aq'^. 

In order to extend Eq. [13] to finite temperatures (T), 
we consider the scaling form of the correlation length, 
^ = VAlgl-TlT/Tcoh], where F[t] is a scaling function 
and Tcoh = I^T^- In this case Eq. [13] can be written as. 
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where T]j{g) represents an upper limit for superconduc- 
tivity around the AFQCP. 

The scaling function F(t) has a well known behavior in 
two limiting cases. First, F{t ^ 0) ex 1— + (t^), such 
that F(0) = 1 and for t < 1, i.e., T <C Tcoh this yields 
a Fermi liquid behavior for the staggered susceptibility 
which in the spin fluctuation theory is used to define the 
correlation length^^. Also, neglecting the effect of dan- 
gerously irrelevant interactions (to be discussed below), 
F{t oo) (X t^ where the exponent x is determined 
by the condition that the dependence of the correlation 
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FIG. 1: (Color online) The line qc^(Tr>) — 1 represents an 
upper limit where superconductivity induced by antiferro- 
magnetic spin fluctuations can occur. It's equation is Td = 
I^I^^C""*^ (see equation I14p and is plotted schemati- 

cally in the figure. It separates region I with ^c«f ^ 1 from 
region (II) of local quantum criticality where fluctuations are 
critical in the time directions but local in space. The latter 
gives rise to repulsive, on-site and nearest neighbors interac- 
tions (see text). 



length on g cancels out. This yields x = — 1/z, such that, 
at the quantum critical trajectory (^ = 0, T ^ 0), the 
correlation length diverges as ^ ex T~^l^ . We can also 
show using this asymptotic behavior of the scaling func- 
tion that the point in the phase diagram at which the 
superconducting temperature can attain its maximum 
value is just above the AFQCP, i.e., at ^ = 0. In fact 
dTD/dg (X l^'l^/^ and vanishes at ^ = 0. 

An interpolation formula for the scaling function which 
gives correct results on both limits (t ^ and oo) is given 
by, 
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where we used the value of the exponents u = 1/2 and 
z = 2. Using this expression for the scaling function, Eq. 
[13] can be written as. 
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This has the form of a dome as shown in Fig. [H This 
equation provides a reasonable interpolation for the lines 
of constant correlation length in the region of the phase 
diagram ^ > 0. For negative {g < 0), there may be 
thermal fluctuations that change the scaling behavior of 
the correlation length, as will be discussed below. The 
physical significance of Td is now quite clear. It pro- 
vides an upper limit to the region where superconduc- 
tivity induced by quantum antiferromagnetic spin fluc- 
tuations can exist. While for T this may provide a 
reasonable estimate of the actual superconducting region 
and its shape, for larger temperatures thermal fluctua- 
tions should reduce the critical temperature to values 
well below Td. 

In the theory of the AFQCP for deff > dc, the quartic 
interaction u is dangerously irrelevant^^ . It determines 
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FIG. 2: (Color online) Line of constant correlation length 
(Qc^ = 1) when thermal or interacting antiferromagnetic spin 
fluctuations are taken into account. 



the shape of the Neel line and changes the value of the 
shift exponent tjj from the expected scaling result, tjj = 
uz, to = z/{d + z - 2) (for the AFQCP in 3d, = 
2/3)^^. The scaling expression for the correlation length 
can be immediately generalized to include the effect of 
the quartic interaction u. It is given by^- 



^ = VA\g{T)rG 
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In this equation g{T) = g- uT^I"^ , such that, g(T^) = 
gives the equation for the Neel line. The scaling function 
G(t) has the following asymptotic behaviors; G(t = 0) = 
1 to reproduce the previous zero temperature results, and 
G(t oo) = The latter guarantees the correct 

behavior close to the critical Neel line giT^) = 0, i.e., ^ = 
Q{T)\g{T))\-^ , with the amplitude Q(T) = V^T^. 
When T Tat, g{T) and the correlation length 
diverges with the thermal correlation length exponent v. 
Assuming-^^, z> = = 1/2, we get at ^ = 0, 
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For g 
^ = 0, 



= 0, this yields a temperature for the dome at 
T^{g = 0) = (Ag^/ii)^^^, to be compared with 
Tnig = 0) = Aq^^ obtained previously. This leads to 
Td = ii(T^)^/^ and since is a small number, we ex- 
pect that in general Td <C T^, such that, non-Gaussian 
fluctuations allow in principle for larger critical super- 
conducting temperatures just above the AFQCP. How- 
ever, in this case the lines of constant ^ satisfying Eq. [13] 
should follow closely the Neel line as in Fig|2l The exper- 
imental results show however that the superconducting 
region has a dome shape^. Then, they seem to imply that 
only purely Gaussian quantum fluctuations are effective 
in pairing the quasi-particles. Also notice that along the 
Neel line, for T 7^ 0, the quartic interaction is a rele- 
vant interaction^^. Then, as the spin fluctuations start 
to interact they apparently loose their efficacy in pairing 
the quasi-particles. At least this is what the experiments 
seem to imply. 
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In some heavy fermion systems as one moves away from 
the AFQCP, for example, applying pressure in the sys- 
tem there is a second superconducting dome^^^^^. This is 
generally attributed to pairing due to charge fluctuations 
associated with a valence transition^^. This second dome 
is larger than that associated with the AFQCP extend- 
ing in a wider region of pressures and temperature. As 
pointed out before in the case of pairing by charge fluctu- 
ations the interaction is attractive when the system is in 
the regime of local quantum criticality. So, we expected 
that in this case, superconductivity can occur in a larger 
region of the phase diagram around the relevant quantum 
critical point, as is in fact observed. 

The energy scale of the magnetic glue that fixes the 
region where superconductivity can exist is given by Aq'^ . 
A is the stiffness of the spin fluctuations and Qc a cut- 
off appropriate for a hydrodynamic description of these 
modes. This quantity plays a role similar to the Debye 
energy in BCS superconductors. 

The region of the phase diagram just above the limiting 
superconducting dome is a state of local quantum critical- 
ity. This state is characterized by a single energy scale, 
the coherence temperature, Tcoh = \d\^^ - It has a resis- 
tivity which scales as p ex (T/TcohY ^^r T <^ Tcoh and 
as p (X {T/Tcoh) for T Tcoh- In the case of anisotropic 
lattices, the crossover to the local regime occurs in stages. 
For a tetragonal system with a spectrum of spin fluctua- 
tions given by, axyQ^ -^CixyQy -^clzqI-, as the system moves 
away from the AFQCP it goes from d = 3 to d = 2 and 
finally to = quantum critical behavior. 



The dome shape of the superconducting region seems 
to indicate that the interaction between the spin fluctu- 
ations acts in detriment of superconductivity. It can not 
be excluded that interacting spin fluctuations can still 
provide a pairing mechanism to produce, for example, 
a pseudo-gap state, but not superconductivity. We have 
shown quite generally that if Gaussian quantum spin fluc- 
tuations give rise to superconductivity, the maximum al- 
lowed Tc can be found just above the AFQCP. Finally, 
our results strongly support the proposal that the second 
superconducting dome observed in some heavy fermions 
systems is due to pairing by charge fluctuations. Since 
even in the local quantum regime, charge fluctuations 
give rise to attractive interactions, superconductivity in 
this case can extend over a wider region of the phase 
diagram. 

Our results are appropriate to describe the system in 
the paramagnetic region. In the long range ordered mag- 
netic phase of the diagram there are new excitations, the 
spin waves, associated with transverse modes. Then in 
this region the present approach does not provide any 
insight. 



Acknowledgments 

I would like to thank Dr. P. Monthoux for invalu- 
able discussions. This work is partially supported by the 
Brazilian agencies CNPq and PAPER J. 



* Electronic address: mucio@if.ufT.br I 

^ Philipp Gegenwart, Qimiao Si, Frank Steglich, Nature 

Physics 4, 186 (2008); Tuson Park, F. Ronning, H. Q. 

Yuan, M. B. Salamon, R. Movshovich, J. L. Sarrao, J. D. 

Thompson, Nature 440, 65 (2006). 
^ M. A. Continentino, G. Japiassu, A. Troper, Phys. Rev. B 

39, 9734 (1989); M. A. Continentino, J. de Physique, I, 693 

(1991); M.A. Continentino, Phys. Rev. B 47, 11587 (1993); 

M. A. Continentino, Phys. Rep. 239, 179 (1994); M.A. 

Continentino, Zeitschrift fiir Physik B 101, 197 (1996). 
^ N. D. Mathur, F. M. Grosche, S. R. Julian, I. R. Walker, 

D. M. Freye, R. K. W. Haselwimmer, G. G. Lonzarich, 

Nature 394, 39 (1998). 
^ P. Monthoux, D. Pines and G. G. Lonzarich, Nature, 450, 

1177 (2007). 

^ P. Monthoux, A. V. Baltsky, D. Pines, Phys. Rev. Lett. 67, 
3448 (1991); Phys. Rev. B 46, 14803 (1992); P. Monthoux 
and G.G. Lonzarich, Phys. Rev. B 59, 14 598 (1999). 

^ P. Monthoux and G.G. Lonzarich, Phys. Rev. B 63, 
054529 (2001). 

P. Monthoux and G. G. Lonzarich, Phys. Rev. B 69, 
064517 (2004). 

^ P. Monthoux and G. G. Lonzarich, Phys. Rev. B 66, 
224504 (2002). 

^ Hisashi Kondo, Torn Moriya, larXiv:0809.4825V 2 [cond- 
mat.str-el], and J. Phys. Soc. Jpn. 78, 013704 (2009); Tet- 
suya Takimoto and Torn Moriya, Phys. Rev. B 66, 134516 



(2002). 

John A. Hertz, Phys. Rev. B14, 1165 (1976). 
""^""^ M. A. Continentino, Quantum Scaling in Many-Body 
Systems^ WorldScientic, Singapore, (2001); M. A. Conti- 
nentino, Phys. Rev. B 47, 11587 (1993). 
T. Moriya, in Spin Fluctuations in Itinerant Electron 
Magnetism ( Springer- Ver lag, Berlin Heidelberg, 1985); T. 
Moriya, T. Takimoto, J. Phys. Soc. Jap. 64, 960 (1995). 
M. A. Continentino, Eur. Phys. J. B 13, 31 (2000); ibid, 
JMMM 272-276, 231 (2004); ibid, Phys. Rev. B 57, 5966 
(1998). 

C. Pinettes and C. Lacroix, J. Phys.-Cond. Matter 6, 10093 

(1994) ; C. Lacroix and C. Pinettes, Physica B 206-207, 11 

(1995) . 

A. J. Millis, Phys. Rev. B 48, 7183 (1993). 
In the original spin fluctuation theory of the QCP^^, the 
coherence temperature is never explicitly introduced. How- 
ever, it is clearly there, with the same significance that it 
was first introduced in Ref.^, as marking the onset of Fermi 
liquid behavior. 

J. Flouquet in W. L Halperin (Ed.), Progress in Low 
Temperature Physics, Vol. 15, p. 139, Elsevier, Amsterdam 
(2005). 

A. T. Holmes, D. Jaccard and K. Miyake, Phys. Rev. B 
69, 024508 (2004). 

A. Kopp and S. Chakravarty, Nature Physics 1, 53 (2005). 



